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Abstract 

In this article we address the problem of quantum tunneling of 
a non-Markovian Brownian particle away from thermal equilibrium. 
We calculate the Kramers escape rate at low temperature (including 
the zero temperature case) in the Smoluchowski limit (strong friction 
regime). Our main findings are: (i) our quantum escape rate is valid 
far from the thermal equilibrium and is non-Markovian, but it becomes 
Markovian as the correlation time vanishes; (ii) at thermodynamic 
equilibrium we obtain a non-Markovian quantum rate that predicts 
a superfluidity phenomenon in the Markovian limit at low and zero 
temperatures. 
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I Introduction: Quantum tunneling 



Let us consider a quantum particle moving in a double-well potential. The 
transition of this particle over the barrier potential from a metastable state 
toward another state is known as quantum tunneling. This phenomenon 
takes place in many areas of physics (e.g., in condensed matter physics), 
chemistry, astronomy, and biology [H [2j [3]. From the experimental stand- 
point quantum tunneling has been investigated on the basis of rate experi- 
ment performed in the following areas [2]: Tunneling in biophysical transport, 
quantum diffusion in solids, chemical conversion processes, tunneling in fer- 
romagnetic materials, electron tunneling in amorphous alloys, nucleation of 
vortices in Hell, escape of electrons from the surface of liquid helium, low- 
temperature Josephson-j unction systems, tunneling of protons on hydrated 
protein powders, and many others. 

In isolated systems quantum tunneling process has been studied making 
use of the Schrodinger equation [I], whereas in open systems, i.e., systems 
immersed in a reservoir undergoing Brownian motion, we cannot follow such 
approach since there is no wavefunction describing a Brownian particle, for 
instance. The nontrivial issue is therefore to find out a manner of elaborating 
a theory of quantum transport so as to clarify the role played by dissipation 
and fluctuation on quantum tunneling. To be complete this theory has to 
take into account not only equilibrium and Markovian properties but also 
non-Markovian and non-equilibrium features [5]. 

The usual Hamiltonian approach for describing open systems is as fol- 
lows. Imagine an environment consisting of a set of harmonic oscillators 
coulpled to the Brownian particle. Once found the Hamiltonian function of 
whole system (particle plus environment), canonical quantization procedure 
could be employed. The Feynman path integral formalism is used to derive 
equations of motion for the reduced density matrix describing the quantum 
motion of the Brownian particle alone. On the basis of this method Caldeira 
and Leggett j6] have quantized the Kramers equation [7j and analyzed quan- 
tum tunneling [5|. More recently, following the Caldeira-Leggett approach 
Ankerhold at al. [9J have derived a quantum Smoluchowski equation and 
explored its physical meaning applying it to chemical reactions, mesoscopic 
physics, and charge transfer in molecules. Meanwhile, other works [TU] have 
revealed some drawbacks underlying such quantum Smoluchowski equation, 
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since it violates the Second Law of Thermodynamics, for instance. From 
our point of view the restlessness over this internal inconsistency leads us to 
search for alternative ways of quantizing the Smoluchowski equation, thereby 
eschewing any ad hoc procedures as achieved in Ref.[10j. 




Another Hamiltonian approach [IT] has been developed and used in [12] 
to derive a quantum Smoluchowski equation and investigate the issue of 
quantum tunneling at zero temperature, in contrast with Ankerhold et al.'s 
survey. Such alternative way does not depend on the path integral formalism, 
but is based on the canonical quantization. 

In recent years, we have put forward a non-Hamiltonian method for quan- 
tizing open systems [131 HU US]- There we start directly with the stochastic 
dynamics (Langevin and Fokker-Planck equations) and quantize it by mak- 
ing use of a Fourier transform carrying the Planck constant. Accordingly, 
our approach is independent of path integral techniques and is not based on 
canonical quantization. 

In order to contribute to a general theory of quantum tunneling within a 
non-Hamiltonian framework we have already started in [15] we organize our 
article as follows: 

I. Introduction: Quantum tunneling 

II. Our generalized Langevin equation 

III. Our non-Markovian Smoluchowski equation 

IV. Quantization of our Smoluchowski equation 

V. Quantum tunneling 

VI. Summary and discussions 

Appendix: Derivation of our Fokker-Planck equation [Eq.(13)] 

II Our generalized Langevin equation 

As a physical model of a stochastic process we consider a particle with mass 
m immersed into an environment. This particle undergoing a Brownian mo- 
tion is characterized by the stochastic position X = X(t) and the stochastic 
momentum P = P(t), while the environment is specified by a random vari- 
able ^ = ^(t). Such physical quantities could be intertwined through the 
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relations 

dX 

X = Q + AQ ; p = ™-ft> i 1 ) 

where AQ = abi(t)^$?(t), t being a parameter, called time, and a a dimen- 
sional constant such that AQ has dimension of length, d/dt denotes a differ- 
ential operator acting upon X, and b\ (t) a time-dependent parameter mea- 
suring the strength of the environment effects upon the particle. We define 
it as being 

h = b 1 {t)= f (*(t'Wt"))dt", (2) 
Jo 

where the mean 

(*(*W)> = J f JW)1>(1?')D X pv(x,p,1>,t)dxd P dil> = 

^(t')^(t")Dy(^,t)dij 

is calculated in terms of the joint probability density function D XP ^(x, p, ip, t) 
or the probability density D^(ip,t). 

One assumes the motion of the Brownian particle moving in an external 
potential V(X) to be described by the stochastic differential equations in 
phase space (X, P), known as Langevin's equations [21 E], 

dP dV 7 „ , T dX P 
dt dX m dt m 

where —jP/ m denotes a (memoryless) frictional force activating the particle 
motion. There $ has the statistical properties 

(*(t')*(t")> = 2D 1/3 5(t" - 1') ■ (*) = 0, (4) 

making the stochastic process Markovian. 5(t"—t') is the Dirac delta function 
and D is a constant - to be determined by the physics of the problem - such 
that bi^ = D 1 / 3 in Eq.(3) has in fact dimension of newton. 

It is important to note that as the environmental parameter &i(i) does 
vanish, the stochastic quantities P and X reduce to the respective determin- 
istic values p = mdq/dt and x = q, provided D XP (x,p) = S(x — q)5(p — p'). 
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Physically, that means that the initially open system becomes isolated from 
its environment and turns out to be described by Newton's equations 

dp dV(x) p dx p , . 

dt dx ^ m ' dt m 

For this reason one says that the Langevin equations (3) are a generalization 
of Newton's equations (5). 

In the literature (2] the non-Markovian character is introduced by means 
of the following statistical properties of ^ 

= PAc) 1/3 e~ (t "~*' )Ac ; (*> = 0, (6) 

where t" > t! and t c is the correlation time between the Brownian particle 
and its environment. One takes into account a memory friction kernel in the 
Langevin equations (3): 

dP dV f l a . .P(t) , T dX P 

(3(t-T)^^dT + b^ ; — = — • (7) 



dt dX Jq m dt m 

Both the frictional kernel (3{t — r) and the fluctuating function are 
coupled by means of the dissipation-fluctuation theorem [16] 



(y(t')y(t")} = K B Tf3(t-r). 

Physically, such a theorem assures that the Brownian particle will always 
attain the thermal equilibrium of the heat bath characterized by Boltzmann's 
constant Kb and the temperature T. As /3(t — r) = 2 , -f5(t — r) and the 
correlation time t c tends to zero, i.e., t c — > 0, the expression (6) reduces to 
(4) while (7) reproduces (3). Thereby, the stochastic dynamics (7), along with 
the statistical properties (6), are called the generalized Langevin equations 
1161. 



In the present paper our purpose is to make another extension of the 
Langevin approach. To begin with, we hold the definition of X in (1) and 
generalize the stochastic momentum P = dX/dt according to 

P = P + AP, (8) 

where AP = — m6 2 (t)\l / (t), with 6 2 (^) defined as 

b 2 = b 2 (t) = f\w))dt'. (9) 
J o 
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Accordingly, the Langevin equations (3) turn out to be written as 

dP dV 7 - , T dX P , T 
— = -—-±P + b 1 y ; — = - + b 2 V, 10 
at dX m dt m 

in phase space (X, P) , with 

(tt(i')*(t")> = (£>Ae) 1/3 e" ( *""*' )Ac ; (*) = (C/^) 1 /3 e -t'A C) 
and 

6 X = (^^(i _ e -t/tc) . fe2 = (Ct c ) 1 /3( 1 _ e -Vtc)_ ( 12 ) 

As the constant C vanishes, we recover from (10) the usual Langevin 
equations (3) as a special case. In short, equations in (10), together with 
(11) and (12), are our generalized Langevin equations. 



Ill Our non-Markovian Smoluchowski equa- 
tion 

Equations (10), (11) and (12) generate the following Fokker-Planck equation 
in phase space (x,p) (for details, see Appendix) 

dJF_ = d{A x T) _ d(A p J r ) + A^(PF_ + A d 2 T + A^d 2 ^ 
dt dx dp 2 dx 2 xp dxdp 2 dp 2 ' 

where 

F = F(x,p,t) = J D xP ^(x,p,ip,t)dip. 

The quantities 

A x = {p/m) + (C 2 /g 1 /3 (e -t/*c _ e -2t/ tc)) 

and 

A, = - hMp+iCD/ty^ie-^ - e- 2t ' tc ) 

dx 

are the drift coefficients, whereas the time-dependent diffusion coefficients 
are given by 

A xx = (C 2 D)^(l-e-^) 2 , 
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A xp ={D*Cy/\l-e-^)\ 

and 

Afp = D(l - e-^f. 
Combining A xx , A xp , and App we notice that they satisfy the relation 

\J A xx App A x p. (14) 

Moreover, on replacing the Maxwell-Boltzmann (MB) distribution 

JT( X p) = 1 e -( P y2m KB T) e -(k X */2 KB T) ^ 

into our Fokker-Planck equation (13) it is too easy to verify that (15) cannot 
become its solution. This means that our stochastic process, described by 
(10-13), holds always away from the thermal equilibrium. That leads us to 
think that the physical meaning of the relation (14), which is a consequence 
of our assumption (\&) ^ in (11), is connected with nonequilibrium char- 
acteristics underlying the environment. In fact, as C = the constraint (14) 
is broken up and our generalized momentum P in Eq.(8) equals to P. Con- 
sequently, Eq.(13) reduces to the non-Markovian Kramers equation in phase 
space (x,p) 



dT p dT d 
dt m dx dp 



— + 1 P)T 

dx m ' 



£)(! _ e -t/*c)2 d 2 

+ 2 (16) 



In the Markovian steady regime characterized by t ^> t c , or formally t c — > 0, 
the MB distribution (15) with p = p turns out to be a solution to (16), thereby 
determining the diffusion coefficient as being equal to A pp = D = 2^kbT. It 
is worth noting that according to both the phase space equations (13) and 
(16) the following fluctuating relation is valid 

AXAP > 0, 



where AX = ^(X 2 ) - (X) 2 , and AP = ^(P 2 ) - (P) 2 . 

On the other hand, inserting F{x,p,t) = f(x,t)S(p) into (13) and taking 
into account the high friction condition 

p dV 
m dx 
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obtained from Newton's equations (5) on neglecting inertial effects (\dp/ dt\ <C 
|7p/m|), we arrive at the non-Markovian Smoluchowski equation in position 
space 

df(x,t) 1 d r „. . A xx d 2 f(xA) . . 

where 

K.(x, t) = ~ + 7 (C 2 /t c ) 1 / 3 (e-^ - e" 2 ^). 
dx 

Replacing (15) into (17) we obtain A xx = Ik^T ji in both stationary and 
Markovian regimes. 

Considering V = (free particle) and C = from our equation (12) we 
derive the non-Markovian Rayleigh equation in p-space [17] 

dg(p,t) id D(l - er^f d 2 

-^T = mdp- [v9M] + 2 W 9{P ^ (18) 

with 

g(p,t) = / J r (x,p,t)dx. 



At thermal equilibrium we find D pp = ik^T as being the diffusion coefficient 
in momentum space after inserting (15) into (18). For both (16) and (18) 
the fluctuating relation 

AXAP = 

is satisfied. 

From the mathematical viewpoint we note that we can derive the Kramers 
equation (16), the Smoluchowski equation (17), and the Rayleigh equation 
(18) as special cases of our equation of motion (13). That physically means 
that all the physics encapsulated into these equations of motion (16), (17), 
and (18) are in principle contained in our Eq.(13). 

In the next sections we wish to survey the quantization of our Smolu- 
chowki equation (17), thereby tackling the problem of quantum tunneling at 
zero temperature in the strong friction regime. 
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IV Quantization of our Smoluchowski equa- 
tion 

We start with the non-Markovian Smoluchowski equation (17) at points X\ 
and x 2 ; we multiply the first equation by f(x 2 ,t) and the second one by 
f(xi,t). We add the resulting equations and obtain the following evolution 
equation for the function £ = ^(xi,x 2 ,t) = f(xi,t)f(x 2 ,t): 

' JC(x 1 ,t) + ^-JC(x 2 ,t))^ + 2A xx ^. (19) 



at 7 \ axi ax2 / 

We now perform the change of variables 

£i + x 2 

x = — - — , n = xi-x 2 . 

We define our quantization process by introducing the following Fourier 
transform [131 M ED 



F(x, p,t) = —j £(x, rj, t)e^ h d v , (20) 

h being the Planck constant h divided by 2tt and p = mdx/dt the physical 
momentum. Inserting (20) into (19) we arrive at our non-Gaussian, non- 
Markovian equation of motion in quantum phase space (x, p; H) 

^ = -h OF+M r +BF)+ m^- mt) p, ( 21 ) 



with 



OF = 2 V 1 d^t)d^F 

^ (s- l)!^) 8 " 1 dx s dp s ~ 1 ' K } 

s=l,3,5,... ' 



oo 



and 



v i dnc^d^F 

r=i?4 r!2r dxdpT ' 

V" 1 dnCi^t) j d s - l F d s F\ 



V(t) = (DC 2 ) 1/3 (1 - e^ /4c ) 2 . (25) 
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Our quantum Smoluchowski equation (21) is valid far from the thermal equi- 
librium and for nonlinear external forces. We emphasize that according to 
our quantization procedure, based on the Fourier transformation (20), our 
Eq.(21) has to obey the Heisenberg fluctuating relation 

AXAP > -, 
~ 2' 

thus restoring the stochastic character of momentum variable in the quantum 
domain. In the classical limit, h — ■> 0, we recover the classical expression 
AXAP > 0. Most specifically, assuming 

F(x,p,t) = f(x,t)5(p) 

our quantum equation (21) leads to the classical Smoluchowski equation (17), 
since 

<9/C 



J OFdp = 2—f(x,t), 
df(x,t) 



i 
I 



AFdp = /C- 



dx 



BFdp=~f(x,t). 



IV. 1 Harmonic oscillator 

Let us consider V = kx 2 /2, k being a constant. Equation (21) turns out to 
be written as 

9F k k r , ...dF V(t)d 2 F k OF _ .p 2 ^ 
W =-F + -[ X+ hit)]- + - -p- - 2V(t/ ¥ F, (26) 

where ^ 

With the initial condition 

F(x, p,t = 0) = y^ e -«P 2 /^ 2 /\ (27) 
ixh 
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a and b being Gaussian parameters, the solution for (26) reads 

1/2 



F ( T „ f) - J_ ( p -A(t)p 2 -[x-c(t)] 2 /4B(t) 



where 



and 



C\ 2/3 fe" 2 */* 



t, 



m = ( f ] < — - - e-^ [> . (33) 



Solution (28) leads to 



that is, 



(28) 



s( * ) = \^ _ 2 ; e + T' (30) 

c(t) = (l-e^/>(t), (31) 

, (t) = 7(c2c) V3{^_ 2tc _^_ + ^L_}, (32) 



AX = v^t), (34) 
1 



AP = 1 , (35) 



V Quantum tunneling 

Now, let t = At be a fixed time interval for observing the Brownian particle 
such that we have a steady solution for (26), that is, dF/dt\ t= A T = 0. During 
At the system is therefore stationary. In this context, we aim at to calculate 
the quantum Kramers escape rate of a Brownian particle over a potential 
barrier in the strong friction regime. 

We consider a Brownian particle moving in a double- well potential V(x). 
The barrier top is located at point Xb, while the two bottom wells is at x a 
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and x c , such that V{x a ) = V(x c ) = 0, x a < x^. The starting point is our 
solution (28) we modify according to 

F(x,p, At) = a( f)(x,p)e- A{AT)p2 - lx -< AT)]2/4BiAT \ a = constant. (37) 

Inserting (37) into the steady version of (26) we derive the ordinary differen- 
tial equation for the function <f>(x,p) 

-2kde=^Te i=p ~ x > (38) 

since 

B<Ar) = % (39) 

A(Ar) = g, (40) 

h(Ar) = -c(Ar). (41) 

From (39) and (40) we determine the parameters of Gaussian function (27) 
as being 

(42) 



and ^_ 

6 = ^{i +2M (l-e^/7)}' (43) 

with 

(- p -Ar/t c f -2Ar/t c >> 

« = ^) = {- 2 '^ + 2k^)}' <44) 
It follows then that ab — 1, whereas from the identity (41) we have the 
following relation among the time scales At (the observation time), t c (the 
correlation time), and t r = j/k (the relaxation time): 



2t r (e-Ar/tc.x)- ^ 
A solution to the differential equation (38) is given by 



-Ar/t 



- 1) (e-^ - 2) 



#0 = J-^ / e (fc/C7)52 ^, A;<0, (46) 



oo 
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wherein we have used the boundary condition 0(£ — > +00) = 1. This result 
confirms the fact that the region around the barrier at Xj,, in which the 
curvature of the potential is negative, is quite relevant in the calculation of 
the diffusion current, as we will see below. Substituting (46), (39), and (40) 
into (37) and expanding the ensuing F(x,p) around Xb we obtain 



F = aC- 2V(Xb)/Vy J c -(V^/nuvih2)pl+(rn^/Vy)(x-x h ~cf f ^{mwl/V^ fa 

V^7 J -00 

(47) 

with kb = -mw k 2 . At point x = Xb we find the following diffusion current 

F(x = x b ,p)^dp = aU f 6 =e -2V( Xb )/v, (48) 

where 

c = c (Ar) = (1 - e m <^h) ^ j r_ e ~^/t c I (4Q) 

Our result (48) is valid provided 

mv%h > P7. (50) 



At the vicinity of x a we cannot use the stationary solution (47), since it 
is only valid for negative curvature k < 0, hence we use the (nonnormalized) 
function 

F(x,p) = ae -(^/^-(k/v l){x -cf 
to finding the number of Brownian particles injected around x a \ 



^=S S F{x ' p)dxdp=aA 



Using v a and the current (48) the non-Markovian quantum Kramers escape 
rate at non-equilibrium regime reads 

• 2fe rrwj 2 c 2 /V-y 

r = h. = ^ fte p -2v(x b )/vj / 51 x 

v a 27r£> 7v /l - {V-i/muolK) 2 
13 



In the Markovian limit At ^> t c , or formally t c — > 0, from (51) we obtain 

T = ^ _ P -2V( Xb )/v, (52) 

2ttP 7v /1 - (V^/muolK) 2 

Non-Markovian properties are therefore responsible for the enhancement of 
the quantum tunneling rate far from the thermal equilibrium. 

We want now to evaluate the diffusion coefficient T> present in (51) or (52) 
using the well-established principles of equilibrium thermodynamics. To this 
end, let us assume that during At our open system has attained a thermal 
equilibrium situation in which is valid the principle of energy equipartition 

that associates the stochastic dynamics of the Brownian particle (e.g., its 
average total energy) to equilibrium thermodynamics underlying the thermal 
environment (Boltzmann's constant Kb and temperature T). 

Replacing (39) and (40) into (34) and (35), respectively, we find 

= |f , (54) 

(P 2 ) = ^ (55) 
that lead to the following equation according to (53) 

7 m7 z 

whose solution is the quantum diffusion coefficient [T5] 

V = + -^(2k b T) 2 - (htt) 2 ; n 2 = k/m, k > 0, (57) 

7 7 

valid for any finite temperature T > HQ/2kb, or 

V = <2 ^- + -y/(2 KB T) 2 + (M 2 ; oo 2 = -k/m, k<0, (58) 

7 7 
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valid for any finite temperature T > 0. At high temperature, KbT ^> HQ, fuv, 
both (57) and (58) lead to the classical diffusion coefficient T> = AkbT/'j. 
On the other hand, at low temperature so that the quantum energy of the 
Brownian particle is equal to the thermal energy of the reservoir, i.e., kbT = 
hVL/2 = tku/2, (57) and (58) lead to V = hVL^ and V = Hlu(1 + \/2)/7, 
respectively. At lower temperature, KbT -C hw, from (58) we obtain D = 
(ftw + 2k_bT)/7 which in turn leads to T> — fku/j at zero temperature. 

Due to the condition (50) the non-Markovian quantum Kramers escape 
rate is thereby constrained to the low-temperature realm (T < f\w b / Vb) 

,2fc iridic 2 /(uj b h+2K B T) 

Y = uj b ne e -2V(x b )/(u b h+2K B T) ^ 

2it{u b K+2K B T)y/\ - [(uj b h+2K B T)/mujlh} 2 

that in the Markovian limit leads to 

e -2V(x b )/(u b h+2K B T) _ ^ 



2ix(uu b h + 2k b T) y/l - {(u b h + 2k B T)/mujlh} 2 
At zero temperature from (59) we obtain the non-Markovian rate 

, , mui b c 2 /h 

r = u b e p -2V( Xb )/u b n (61) 

2ny/l - {l/mw b ) 2 
whereas in the Markovian case we derive 

r - ^ p -2v(x b )/u b n 



e -2V(x b )/u b n ^Q2) 



2n^l - (l/muj b ) 2 
that is the same result obtained in our previous work [15] . 

We wish to point out that our Markovian results (60) and (62) are inde- 
pendent of the friction constant 7 at thermal equilibrium. This means that 
over the barrier at low temperatures (including zero temperature) Brownian 
particles may decay into a metastable well around x c at high temperatures. 
This may take place without violating the Second Law of Thermodynamics 
since in the quantum domain at low temperatures the particles could flux 
overcoming any dissipative mechanism (superfluidity phenomenon!). 

So far we have considered the thermal environment as having a general 
nature. Supposing then that the thermal reservoir consists of a set of many 
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harmonic oscillators we can generalize the equipartition theorem of energy 
(53) taking into account the quantum nature of the heat bath: 

(E) = + -(X 2 ) = (hu/2) coth(hv/2K B T). (63) 
2m 2 

[u = (k/m) 1/2 for k > 0, and v = (-k/m) 1/2 for k < 0]. It follows then the 
equation V 2 — A(E)V/'y + kh 2 jm^ 2 = that yields the solution 



^ fiU fiQ Vth I . 9 fin ^ 9 , . , „ . . 

V = — coth + — Wcoth 2 1; tt 2 = km, k > 0, (64) 

7 2k b T 7 V 2k b T 



or 



^ fiu , fiw uoh . 9 ftw 9 , . , , . 

V= — coth- — - + — \/coth 2 - — - + 1; iu 2 = -km, k<0. 65 
7 2n B T 7 V 2k b T 

At high temperature, k b T ^> Ml, Tiu;, both (64) and (65) lead to T> — 
4K B T/ / ~f. At low temperature, k b T <C Ml, fiio, (including zero temperature 
T = 0), (64) and (65) lead to V = Ml/'y and v = Ml + V2)/7, respectively 
Hence we obtain the non-Markovian quantum Kramers escape rate (52) as 
being given by the expression 

r = ^ __ e -2V(x b )/hw b (l+V2) _ ^ 



2tt(1 + V2)^/l- [(l + V2)/mw b } 2 



In the Markovian limit we obtain 

r = Ub p -2V(x b )/hu, b (l+y/2) ^ 

2tt(1 + y/2)yjl- [{l + V2)/mutf 

that exhibits the influence of the zero-point energy of the thermal reservoir 
upon the quantum tunneling rate of our Brownian particle at zero temper- 
ature. Eq.(67) is also independent of frictional constant 7, thus leading to 
the superfluidity phenomenon. 



VI Summary and discussions 

In this paper we have addressed the issue of quantum tunneling at low (in- 
cluding zero temperature) for non-Markovian open systems away from the 
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thermal equilibrium. From our quantum Smoluchowski equation (21) we 
have derived the non-equilibrium, non-Markovian quantum escape rate (51) 
that does depend on both the friction and diffusion coefficients. 

We have thus provided an alternative method of quantizing the Smolu- 
chowski equation taking into account non-Markovian and non-equilibrium 
effects. This result, coming from our non-Hamiltonian account, is novel and 
is not present in others approaches [HI E21 ESI ED] ■ 

At low temperature and in the Markovian regime our quantum Kramers 
escape rates (60) and (67) are independent of damping constant 7. That is, 
dissipation has no effect on the quantum tunneling, thus giving rise to su- 
perfluidity. This result has been misunderstood by Ankerhold et al. [19] that 
state: "This finding (...) contracts existing theoretical results verified by ex- 
perimental data" . This statement is incorrect since Ao et al. [21] have arrived 
at our same conclusion by studying Landau-Zener tunneling in a dissipative 
environment: "no effect of dissipation is present at zero temperature". 

We hope our present work could foster development of experimental re- 
searches on quantum tunneling in overdamped open systems in order that 
we can compare our theoretical predictions with experimental data. 

In a forthcoming paper [22] we pursue our non-Hamiltonian approach 
in surveying quantum tunneling effects from the quantization of our non- 
equilibrium, non-Markovian Fokker-Planck equation (13) and of the Rayleigh 
equation (18), as well as quantum tunneling in the small friction regime 
(energy diffusion). 

To conclude, we want to point out that non-linear effects due to the 
external potential can be also analyzed through our quantum Smoluchowski 
equation (21) 

dF_ _ 1 / 2 dF_ + 1 d 3 F + ^ F _ 2x dF\ + V(t) d 2 F _ 2V ^P 2 F 
dt 7 \ dx 4 dxdp 2 dp J 2 dx 2 h 2 

(68) 

for the cubic potential V = x 3 /3, for instance. 
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Appendix: Derivation of our Fokker— Planck 
equation [Eq.(13)] 

In this appendix we want to show in somewhat details how we could explicitly 
construct the Fokker-Planck equation (13) from the system of stochastic 
differential equations 



— = -—-±P + b 1 * ; =- + b 2 m. (69) 

at dX m at m 

Equations (69) yield the results 

AP = ~ ( Jx + m P ) At + MO W)>^' (70) 

and 

ax =iu - ( % + ip) ^s! + 1 r Ai r W)mmtis + 

m \ dX m J m m J t J t 

t+At 

b 2 {t')^{t'))dt'.{n) 



Using AP = P(t + At) - P(t) and AX = X(t + At) - X(t) we calculate the 
following quantities 



A — lim — ^ 

Ap ~A^o At ~ dX 



-l p+ ^hST mw!)) ^ (72) 



AM = lim = -2 ( + lp) lim F * l J 1 (0(*(0)*' + 

At~*o At \ dX m J At^o J t 



t+At r t+At 

l\U fj.H\ /VW-lAVTr/j/'W J 4-1 J 4-1 1 



Hm Q ^ / / hdt')bdf){Wmt"))Mdi f '.{TZ) 



t Jt 



A x = lim = - + - lim / / bi(t'){V(t'))dt'ds + 

At^o At m mAt->oJ t J t iV A v ;/ 

& 2 (*')W)X (74) 



At^O 
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A xx = lim ((A y = h + h + h + h (75) 

At-t-0 At 



with 

r t+At rs 



np rt+ixt rs 

h = — / h(lf)(W))dlfds, (76) 

m Jt Jt 

2p rt+At 

h = — b 2 (t')(*(t'))dt', (77) 

m Jt 

-i rt+At rt+At rr rs 

7 3 = -^lim / / / / b 1 (t r )b 1 (t")(^(t')^(t"))dt'dt"drds, (78) 
m At^o J t J t J t J t 



r> rt+At rt+At rs 

7 4 = — Urn / / / b 1 (t')b 2 (t")(y(t')y(t"))dt'dt"ds, (79) 
"' A; ".// .// A 

/•t+At /-t+At 

/ 5 = / / ^(OMOWO^"))^'', (so) 
jt J* 

and 

Axp = 21™ ^ = 6 + 6 + 6 + & + 6, (81) 

where 

p /-t+At 

6 = -/ (82) 

i / jy \ rt+At rs 

i i rt+At rt+At rs 

e 3 = -lim — / / / 6 1 (06i(O(*(0*(O>*'* ff ^ (84) 
m At-»o zit 7* Jt Jt 

p rt+At 

& = / b 2 (t'){V(t'))dt' (85) 

m J t 

i rt+At rt+At 

6 = iSo At i i Ht'Mf)(mW'))dt'dt". (86) 



After using our definitions 
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(V(t')V(t")) = (^Ae) 1/3 e^ ( '"~*' )Ac ; (*) = (C/^) 1/3 e-'' /tc , (87) 

and t 

fe, = /" = (^^(l - e-^), (88) 

Jo 

6 2 = [\y(t'))dt' = (Ct c y/ 3 (1 - e-'/' c ). (89) 
Jo 

into (72-75) and (81) we obtain our Fokker-Planck equation (13) with the 
coefficients 

!_j ( e -*/tc_ e -2t/t e)> (90 ) 

A p - = -jfex - ( + ( ^y /3 (e~* Ac - e" 2 * Ac ), (91) 

A xx = (C7 2 £)) 1 / 3 (l- e -*/tc)2 > (92) 

A xp =(D 2 C) 1 ' 3 (l-e- t ' t °) 2 , (93) 

^ = D(l - e-^f. (94) 
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